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Like early-type galaxies, also nearby galaxy clusters define a Fundamental Plane, a luminosity- 

^ . radius, and a luminosity-velocity dispersion relations, whose physical origin is still unclear. By 

means of high resolution N-body simulations of massive dark matter halos in a ACDM cosmology, 

we find that scaling relations similar to those observed for galaxy clusters are already defined by 

their dark matter hosts. The slopes however are not the same, and among the various possibilities 

in principle able to bring the simulated and the observed scaling relations in mutual agreement, 

' we show that the preferred solution is a luminosity dependent mass-to-light ratio {M/L oc 

. that well corresponds to what inferred observationally. We then show that at galactic scales there 

is a confiict between the cosmological predictions of structure formation, the observed trend of 

^ • the mass-to-light ratio in ellipticals, and the slope of their luminosity-velocity dispersion relation 

(that significantly differs from the analogous one followed by clusters). The conclusion is that 

the scaling laws of elliptical galaxies might be the combined result of the cosmological collapse 

of density fluctuations at the epoch when galactic scales became non-linear, plus important 

03 . modifications afterward due to early-time dissipative merging. Finally, we briefly discuss the 

possible evolution of the cluster scaling relations with redshift. 

Introduction also useful to get insights on the structure and, 

possibly, on the formation and evolutionary pro- 

■ Early-type galaxies are known to follow well de- messes of galaxies and galaxy clusters. 

^ , fined scaling relations involving their main obser- ttt n i ^ i i- i ,• ,i , n .i 

■ -r rr • Well dcfancd scalmg relations, that recall the 

vational properties, i.e, the luminosity L, effective , , • i i , i .i i ■ r 

,. „ , , ■ 1- • • • observed ones, are indeed expected on the basis of 

radius Rc, and velocity dispersion a: m partic- . , , j i r , i r , • r , , 

, 1 , x^, Ti n r the simplest model for the formation of structures 

ular we recall here the i<aber- Jackson hereafter . i- tt • i 

1 DTI ,1 , in an expanding Universe, namely the gravita- 

FJ; Faber & Jackson 1976), the Kormendy (Kor- ^. , a ^ ^- ■ 

, 1 1 1 1 Til /I tional collapse of density fluctuations m an other- 

mendy 1977), and the lundamental Plane here- . , i- , -i ,• r n- • i i i 

_ "L^ , . „ . 1 1 wise homogeneous distribution of collisionless dark 

after IP; Diorgovski & Davis 1987; Dressier et al. /T^nT^ Tr^-^i i-w ii. ii 

. T,r- 1 ■ 1 1- ■ • <• matter (DM). In fact, the spherical top-hat model 

1987) relations. Withm the limitations of a poorer c r-i mnn\ j- ^ ii ^ i 

. . r J (Gunn & Gott 1972) predicts that, at any given 
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statistics, analogous relations have been found to 



epoch, all the existing DM halos have just col- 



hold also for galaxy clusters (Schaeffer et al. 1993, , , 9/^/1 

1 n r.r>o A i • . i i r^r^r, 1 n A lapscd aud viriahzcd, i.c. , = Tvir CTvir / ^-^ ( wffere 
hereafter S93; Adami et al. 1998a, hereafter A98; 7^1/^7.^2 2 nrn/,r tt j >n u 

„ „ ' ,. , ^- • , • V ''vir = -U GM\ ct4 = 2T M, U and T be- 

lor 3-parameters scaling relations involving A-ray . ; .■ , V i- ■• • r 11 

, , , . .. „ m , 1 mg the potential and kinetic energies of a halo 

observables, see Annis 1994; i'ujita & lakahara 



of mass M, respectively). In addition, all the 

^ 1 n • 1 1 ■ • 1 • halos arc charactcrlzed by a constant mean den- 

Gorman 1999 . Besides their potential importance ., . , ,, -.-ii •; r;i tt- 

. . sity Pa J given by the critical density of the Uni- 



1999a; Fritsch & Burchert 1999; Miller, Melott & 
Gorman 1999). Besides their potential importance 
as distance indicators, these scaling relations are 



verse at that redshift times a factor A depending 
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on z and on the given cosmology (sec, e.g., Pee- 
bles 1980; Eke, Cole & Frank 1996). For simplic- 
ity, we call rA the radius of the sphere contain- 
ing such a mean density, so that M oc r\. In 
general, fyir 7^ ^a, but, if for a family of den- 
sity distributions r^ir/r/;^ ~ const, then the virial 
theorem can be rewritten as M oc rAcr^ir, and, 
together with M oc r\, it brings to M oc cr^jj,, 
thus providing three relations that closely resem- 
ble the observed ones. Note that these expec- 
tations involve the global three-dimensional prop- 
erties of DM halos, while the quantities entering 
the observed scaling relations are projected on the 
plane of the sky. However, if DM halos are struc- 
turally homologous^ systems, as found in cosmo- 
logical simulations (Navarro, Frcnk & White 1997, 
hereafter NFW), and are characterized by simi- 
lar velocity dispersion profiles (e.g.. Cole & Lacey 
1996), their projected properties are also expected 
to follow well defined scaling relations (with some 
scatter due to departures from perfect homology 
and sphericity). 

Of course, the simple considerations above are 
not sufficient to account for the observed scaling 
relations of galaxy clusters, at least for two rea- 
sons. The first is that a given potential well (as 
the one associated with the cluster DM distribu- 
tion) can be filled, in principle, by very different 
distributions of "tracers" (such as the galaxies in 
the clusters, from which the scaling relations are 
derived). This means that the very existence of the 
cluster FP implies a remarkable regularity in their 
formation processes: galaxies must have formed or 
"fallen" in all clusters in a similar way. The sec- 
ond reason is that any trend of the cluster mass- 
to-light ratio (necessary to transform masses, in- 
volved in the theoretical relations, into luminosi- 
ties^, entering the observed ones) must be taken 
into account for a proper interpretation of the ob- 
served scaling relations. 

A distinct but strongly related question about 
the origin and the meaning of the scaling laws 
naturally arises when applying the predictions of 

^Strictly speaJsing, DM halos axe only weakly homologous 
systems, since the low mass halos are systematically more 
concentrated than the more massive ones (for a definition 

of weak homology, sec Bcrtin, Ciotti & del Principe 2002). 
^We recall here that the luminosity of a cluster refers to the 
sum of the luminosities of all its constituent galaxies, i.e., 
a sum over the cluster luminosity function. 



the cosmology also at galactic scales: in fact, 
while scale-invariant relations are predicted, dif- 
ferent slopes of the FJ relation are observed for 
galaxies and for galaxy clusters (sec Section 2, 
and Girardi et al. 2000). This suggests that dif- 
ferent processes have been at work in setting or 
modifying the correlations at the two mass scales. 
The theoretical implications of the scaling laws 
for elliptical galaxies (Es) have been intensively 
explored (sec, e.g., Bertin et al. 2002 and refer- 
ences therein), and several works have been de- 
voted to their study within the framework of the 
dissipationlcss merging scenario in Newtonian dy- 
namics (e.g., Capelato, de Carvalho & Carlberg 
1995; Nipoti, Londrillo & Ciotti 2003a, 2003b, 
hereafter NLCOSab; Evstigneeva, Rcshctnikov & 
Sotnikova 2002; Dantas et al. 2003; Gonzalez- 
Garcia & van Albada 2003). In particular, if the 
initial total energy of the system is non-negative, 
the merger products are found to follow the ob- 
served edge-on FP (with the important exception 
of merging dominated by accretion of small galax- 
ies), but they badly fail at reproducing the FJ 
and the Kormendy relations, in accordance with 
elementary predictions based on energy conserva- 
tion and on the virial theorem (NLC03ab). It has 
also been shown (Ciotti & van Albada 2001) that 
gas free mergers cannot account at the same time 
for the FP and for the Mbh-c relation (that link 
the black hole mass and the stellar velocity disper- 
sion of the host spheroid, Fcrrarcsc & Mcrrit 2000; 
Gebhardt et al. 2000); the importance of gas dissi- 
pation in the formation of elliptical galaxies is also 
apparent from the observed color magnitude and 
Mg2-(T relations (e.g., Saglia et al. 2000; Bernardi 
et al. 2003a). Much less effort has been devoted 
to the theoretical study of the FP of galaxy clus- 
ters: beside a work on the effects of dissipation- 
less merging (Pentericci, Ciotti & Renzini 1996), 
the other few theoretical studies mainly focus on 
the possible relation between the FP properties 
and the cluster age, the number of substructures, 
and the underlying cosmology (Fujita & Taka- 
hara 1999b; Beisbart, Valdarnini & Buchert 2001). 
Prom the comparison between the FP of galaxies 
and of galaxy clusters in the /c-spacc, Burstcin et 
al. (1997) derived support to the idea that, at 
variance with groups and clusters, gas dissipation 
must have had an important role on the formation 
and evolution of galaxies. 
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In order to get a more complete view of the 
problems depicted above, we use high-resolution 
N-body simulations to study the scaling relations 
of very massive DM halos, that are thought to host 
the present day galaxy clusters. The aim is to ver- 
ify whether these relations are similar or not to the 
observed ones, and determine the assumptions re- 
quired to make them in mutual agreement. The 
results thus obtained and the empirical evidence of 
different slopes of the scaling relations at galactic 
and cluster scales are then discussed. The paper 
is organized as follows: in Section 2 we present the 
scaling relations observed for nearby galaxy clus- 
ters and those followed by early-type galaxies. In 
Section 3 we describe the high-resolution resimu- 
lation technique employed to build the sample of 
massive DM halos used for our analysis, and we 
derive their scaling relations. In Sections 4 and 
5 we disciiss under which hypothesis these scaling 
relations can be translated into the observed ones, 
focusing in particular on the astrophysical impli- 
cations of the differences at galactic and cluster 
scales. The main results are summarized and dis- 
cussed in Section 6. 

2. Scaling relations of clusters and early- 
type galaxies 

Prom the observational point of view only two 
works (namely, S93 and A98) report on the FP of 
galaxy clusters, i.e., a relation among their opti- 
cal luminosity, scale radius, and velocity disper- 
sion. The two groups agree about the existence 
of a tight and well defined FP, even if quantita- 
tive differences in the numerical values of its co- 
efficients are found, that can be traced back to 
the different choice of the radial variable. In fact, 
A98 use 4 density profiles (King, Hubble, NFW, 
and de Vaucouleurs) to describe the distribution of 
cluster galaxies, concluding that the best fit is pro- 
vided by the first two models: consistently, they 
adopt the clusters core radii to obtain the FP. S93 
instead use in their analysis the cluster half-light 
projected radii (the standard choice in the vast 
majority of the FP studies). For this reason we 
choose to compare our results to those presented 
by S93: we note however that the FP coefficients 
derived by A98 when using the projected half-light 
radii derived from the de Vaucouleurs model agree 
within the errors with those reported by S93. 



For their sample of 16 galaxy clusters at z < 
0.2, S93 used the photometric parameters L and 
Re in the V band (quoted by West, Oemler & 
Dckcl, 1989), and the velocity dispersion a (given 
by Struble & Rood, 1991), and they derived not 
only the FP, but also a FJ-like and a Kormendy- 
likc relations. However, instead of reporting the 
S93 scaling laws, that have been obtained by 
means of least-square fits to the data, here we 
re-derive them by minimizing the distance of the 
residuals perpendicular to a straight line (for the 
FJ and the Kormendy relations) or to a plane (for 
the FP). Results arc anyway consistent with those 
of S93. The FJ and Kormendy relations we obtain 
are: 

LCX^2.18±0.52^ (1) 

in good agreement also with the results of Girardi 
et al. (2000), and 

L oc i?i-55±o.i9^ (2) 

where L is given in 10^^ Lq/H^, a in 1000 km/s, i?e 
in Mpc/h {Ho = 100 /i km s"^ Mpc-\ and h = 1), 
and the errors on the exponents take into account 
also the observational uncertainties. As can be 
seen from Fig.l, where equations (1) and (2) are 
plotted together with the S93 data, the two rela- 
tions above describe real scalings among the clus- 
ter properties, even if their scatter is quite large: 
without taking into account the observational er- 
rors, the rms dispersion of the data around the 
best- fit lines is 0.19 in both cases. Similarly to 
what happens for galaxies, a considerable improve- 
ment is achieved when combining all the three 
observables together in a FP relation, that we 
have derived by performing a Principal Compo- 
nent Analysis (PCA; see e.g. Murtagh & Heck 
1987) on the data sample, thus obtaining the new 
orthogonal variables p, defined by: 

Pi = Q!jlogi?c + Alog-^ + 7»log(7, 1 = 1,2,3. (3) 

The numerical values of the coefficients A, and 
7i are listed in Table 1, while the resulting distri- 
bution of the observed clusters in the (pi,P3) and 
(Pi)P2) spaces are shown in Fig. 2, the former pro- 
viding an exact edge-on view of the FP and mak- 
ing apparent its small thickness. A more intuitive 
representation of a (nearly) edge-on view of the 
FP can be obtained by solving equation (3) for 
i = 3 and using ps ~ const: 

^^^0.9±0.15^1.31±0.22 
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This relation is shown in Fig. 3, and is character- 
ized by an rms dispersion of ~ 0.07. 

How do cluster scaling relations compare with 
the analogous relations followed by early-type 
galaxies? For what concerns the FP, the agree- 
ment with that of galaxies is remarkable. For 
example, the FP of elliptical galaxies in the B 
band is given by L oc cr~^ '^ (e.g., Dressier 

et al. 1987; J0rgensen, Franx & Kjaergaard 1996; 
Scodeggio et al. 1998; Bernardi et al. 2003c). As 
well known, the exponents of the galaxy FP de- 
pend on the adopted photometric band (see, e.g., 
Scodeggio et al. 1998; Treu 2001), but the dif- 
ferences from the values reported above become 
significant only when using the K-band (Pahre, 
Djorgovski & de Carvallio 1998). Thus, the agree- 
ment between the cluster and the galaxy FP can 
be regarded as robust. The situation is similar for 
the Kormendy relation: in fact, L oc 
has been reported for ellipticals in the B band 
(Davies et al. 1983; see also Schade, Barrientos & 
Lopez-Cruz 1997; Bernardi et al. 2003b). How- 
ever, this agreement should be regarded as less 
robust than that of the FP, because largely dif- 
ferent values of the Kormendy slope for various 
selections of the data sample have been reported 
in the case of galaxies (Ziegler et al. 1999). In 
any case, the FJ relation of galaxies is different 
from that the clusters FJ: in fact, L oc ct^ for 
(luminous) Es (Faber & Jackson 1976; Forbes & 
Ponman 1999; Bernardi et al. 2003b), while for 
clusters the reported slope is around 2 (a value of 
~ 1.58 is also found by A98). The additional facts 
that the slope of the galaxy FJ seems to drop be- 
low 3 for ellipticals with a < 170 km/s (see, e.g., 
Davies et al. 1983), and that the exponent of the 
Mbh-(7 relation (Ferrarese & Merrit 2000; Geb- 
hardt et al. 2000) is remarkably similar to that of 
the FJ will be discussed in Section 5. 

3. DM halos scaling relations 
3.1. The simulations 

To investigate whether the dark matter hosts 
of galaxy clusters, as obtained by numerical sim- 
ulations, do define scaling relations similar to the 
observed ones, a large enough sample of very mas- 
sive DM halos is needed. Therefore, we employed 
dissipationless simulation with 512^ particles of 
6.86 X W^^Mq/H mass each, where the volume 



of the Universe is sufficiently large for this pur- 
pose: the box side is 479/i~^Mpc comoving, with 
h = 0.7 (see Yoshida, Sheth & Diaferio 2001). The 
adopted cosmological model is a ACDM Universe 
with = 0.3, Qa = 0.7 (e.g., Ostriker & Stein- 
hardt 1995), spectral shape F = 0.21, and normal- 
ization to the local cluster abundance, (Ts = 0.9. 
From this simulation, we have randomly selected 
a sub-sample of 13 halos at z = 0, with masses be- 
tween 10^^Mrs,/h and 2.3 x lO^^MQ/h. They span 
a variety of shapes, from nearly round to more 
elongated. The richness of their environment also 
changes from case to case, with the less isolated 
halos usually surrounded by pronounced filamen- 
tary structures, containing massive neighbors (up 
to 20% of the selected halo in mass). 

Given the mass resolution of the simulation, less 
than 1500 particles compose a halo of IO^^Mq/Zi 
and, due to discreteness effects, its properties 
defining the FP relation cannot be accurately de- 
termined. We have therefore resimulated at higher 
resolution the halos in our sample by means of 
the technique introduced in Tormen, Bouchet & 
White (1997): here we recall only its relevant as- 
pects (for more details see Lanzoni, Cappi & Ciotti 
2003). The first step is to select in a given cosmo- 
logical simulation the halo one wants to "zoom 
in". Then, the region defined by all the parti- 
cles composing the selected halo and its immediate 
surroundings is detected in the initial conditions of 
the parent simulation, and the number of particles 
within it is increased by the factor needed to attain 
the suited mass resolution. Such a region is there- 
fore called "the high resolution region" (HRR). 
Since the mean inter-particle separation within the 
HRR is smaller than in the parent simulation, the 
corresponding high-frequency modes of the pri- 
mordial fluctuation spectrum are added to those 
on larger scales originally used in the parent sim- 
ulation, and the overall displacement field is also 
modified consequently. At the same time, the dis- 
tribution of surrounding particles is smoothed by 
means of a spherical grid, whose spacing increases 
with the distance from the center: in such a way, 
the original surrounding particles are replaced by 
a smaller number of macroparticles, whose mass 
grows with the distance from the HRR. Thanks 
to this method, even if the number of particles in 
the HRR is increased, the total number of parti- 
cles to be evolved in the simulation remains small 
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Table 1 

pca parameters for the observed clusters. 





Oil 




7i 


Pi 


-0.67 


-0.51 


-0.94 


P2 


0.88 


-0.005 


-1.22 


P3 


0.55 


-0.61 


0.80 



Note. — pi is defined in 
equation(3). 



enough to require reasonable computational costs, 
while the tidal field that the overall particle dis- 
tribution exerts on the HRR remains very close to 
the original one. For the new initial configuration 
thus produced, vacuum boundary conditions are 
adopted, i.e., we assume a vanishing density fluc- 
tuation field outside the spherical distribution of 
particles with diameter equal to the original box 
size L. A new N-body simulation is then nm start- 
ing from these new initial conditions, and allows 
to re-obtain the selected halo at the required res- 
olution. 

We have applied this technique to the selected 
13 massive DM halos. For 8 of them the reso- 
lution has been increased by a factor ^ 33, by 
means of high-resolution particles of mass ~ 2.07 x 
IO^Mq /h each, while a further increase of a factor 
of 2 has been adopted for the 5 intermediate mass 
halos (the particle mass is 10^Mq//i in this case). 
The gravitational softening used for the high- 
resolution region is e = 5kpc//i (roughly Plum- 
mer equivalent), corresponding to about 0.2% and 
0.5% of the virial radius of the most and least mas- 
sive halos, respectively. This scale length repre- 
sents the spatial resolution of the resimulations, to 
be compared with that of 30 kpc//i of the original 
one. To nm the resimulations, the parallel dissi- 
pationless tree-code GADGET (Springel, Yoshida 
& White 2001) has been used. At ^ = the DM 
halos have been selected by means of a spherical 
overdensity criterium (Lacey & Cole 1994; Tor- 
men et al. 1997), i.e., they are defined as spheres 
centered on maximiim density peaks in the parti- 
cle distribution, and with mean density equal to 



the virial density pA predicted by the spherical 
top-hat model for the adopted ACDM cosmology 
(pa — 97 pcrit at z = 0). The corresponding mass, 
linear scales and velocity dispersions are listed in 
Table 2. 

3.2. Scaling relations for the DM halos 

According to the selection method described in 
the previous section, the DM halos arc all charac- 
terized by the same pA and should also be nearly 
virialized systems. The first condition (and thus 
the relation M oc r^) is satisfied by construction 
while the second can be easily verified. In fact, we 
find that all the selected halos follow the relations 
M oc rA T^ir with a rms scatter of 0.03 only, and 
M (X a^ji, with rms ~ 0.05. In Table 2 we list the 
ratio rvir/^A for each halo. 

However, projected quantities are involved in 
the observations and the first step of our anal- 
ysis is the determination of which (if any) scal- 
ing relations are satisfied by the DM halos when 
projected. We have therefore constructed the 
projected radial profiles of the selected halos by 
counting the DM particles within concentric shells 
around the center of mass for three arbitrary or- 
thogonal directions {x, y, and z), and we defined 
Rh as the projected radius of the circle containing 
half of the total number of particles. Then, the ve- 
locity dispersion ah has been computed from the 
line-of-sight (barycentric) velocity of all the parti- 
cles within i?h- Since the DM halos, as well as real 
clusters, are not spherical, such a procedure gives 
different values of Rh and cth for the three line- 
of-sights (the maximum variations however never 
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Table 2 

Properties of the DM halos at 2; = 0. 



Name 


M 




rvir/rA 


tTvir 
















IO^^Mq/Zi 


Mpc/h 




km/s 


Mpc/h 


Mpc/h 


Mpc/h 


km/s 


km/s 


km/s 


1542 


0.82 


0.99 


0.77 


760 


0.22 


0.24 


0.21 


471 


470 


511 


3344 


1.09 


0.99 


0.83 


815 


0.30 


0.24 


0.29 


459 


554 


480 


914 


1.45 


1.09 


0.75 


967 


0.27 


0.28 


0.28 


616 


605 


541 


4478 


2.92 


1.37 


1.02 


1075 


0.54 


0.52 


0.45 


556 


656 


689 


1777 


3.83 


1.50 


0.94 


1177 


0.40 


0.54 


0.46 


783 


638 


742 


564 


4.91 


1.63 


1.01 


1271 


0.65 


0.59 


0.47 


696 


878 


766 


689 


6.08 


1.75 


1.01 


1354 


0.56 


0.65 


0.69 


837 


749 


761 


245 


6.50 


1.79 


1.03 


1332 


0.51 


0.70 


0.76 


855 


799 


764 


51 


10.78 


2.12 


0.84 


1794 


0.62 


0.70 


0.52 


1114 


1018 


1241 


696 


11.37 


2.16 


0.98 


1669 


0.74 


0.72 


0.74 


1062 


1048 


959 


72 


11.77 


2.18 


0.90 


1777 


0.70 


0.71 


0.55 


1105 


993 


1207 


1 


13.99 


2.31 


0.88 


1870 


0.66 


0.72 


0.69 


1187 


1136 


1189 


8 


23.42 


2.75 


0.92 


2262 


0.79 


0.93 


0.99 


1459 


1395 


1319 



Note. — M is the DM halo mass, and rvir arc the truncation and virial radius, respectively, (Tvir is the virial 
velocity dispersion; the projected half-mass radii and the corresponding line-of-sight velocity dispersions are given for 
three orthogonal directions. Note that Rh/i^A is in the range ~ 0.2—0.4. 



exceed 33% and 21% for the two quantities, re- 
spectively; sec Table 2) and wc decided to build 
our data sample by considering all the three pro- 
jections for each halo. With the projected proper- 
ties i?h and (Jh now available, we have determined 
the best fit relations between M and i?h, and be- 
tween M and ah by minimizing the distance of 
the residuals perpendicular to a straight line, and 
thus obtaining the DM analogues of the observed 
FJ and Kormendy relations: 

M oc a^-02±o.i5^ (5) 

and 

M a i?2-36±0.14^ (g) 

with rms ~ 0.12 and 0.15, respectively. With 
a PC A of these data we determined the relation 
analogous to equation (4): 

M oc i?l.l±0.05^1.73±0.04^ (7) 

with rms = 0.04. Compared to those among the 
virial properties, these relations have larger scat- 
ters, as expected. The FJ and FP have slopes 
similar to those obtained for the virial quanti- 
ties, while the M-Rh relation appears to be sig- 
nificantly flatter, as a consequence of the different 
density concentration of low and high mass ha- 
los. Note that the NFW concentration parame- 
ter varies only weakly in the mass range covered 



by the halos in our sample (e.g., Eke, Navarro & 
Steinmctz 2001; Bullock ct al. 2001). This is con- 
sistent with our results: in fact, from equation (6) 
and M oc r|, we find that RiJrA oc M~"-0^. 

We stress that while scaling relations between 
M, ta (or Tvir) and tTvir were expected, a tight 
correlation between projected properties is a less 
trivial result: in fact, structural and dynamical 
non-homology can, in principle, produce signifi- 
cantly different effective radii and projected ve- 
locity dispersion profiles for systems characterized 
by identical M, ryir, and (Tvir- It is also known 
that weak homology, coupled with the virial theo- 
rem, does indeed produce well defined scaling laws 
(see, e.g., Bertin et al. 2002). Therefore, the scal- 
ing relations presented in equations (5) (7) are a 
first interesting result of our study. The differ- 
ence between the values of the exponents appear- 
ing in equations (5)-(7) and those in the virial 
relations is the direct evidence of weak homology 
of the halos: in fact, strong homology would re- 
sult in exactly the same exponents, while a strong 
non-homology would disrupt any correlation. We 
checked this further, and we found that the DM 
halos in our sample still show well defined scal- 
ing relations (although with different exponents) 
when considering the projected radius encircling 
any fixed fraction of the total mass, and the cor- 
responding projected velocity dispersion within it 



6 



(sec Tabic 3). Note that this finding is in agree- 
ment with the results already pointed out by sev- 
eral groups, namely the fact that DM halos ob- 
tained with numerical N-body simulations in stan- 
dard cosmologies are characterized by significant 
structural and dynamical weak homology (e.g., 
Cole & Lacey 1996; NFW; Subramanian, Cen & 
Ostriker 2000, and references therein). 

4. Prom the simulated to the observed 
scaling relations 

Comparison of equations (5), (6), (7) and (1), 
(2), (4), reveals that the FJ, Kormendy and FP re- 
lations of simulated DM halos are characterized by 
different slopes with respect to those derived ob- 
servationally. What kind of regular and system- 
atic trend with cluster mass of the galaxy prop- 
erties and distribution are implied by these dif- 
ferences? In order to answer this question, we 
define the dimensionless quantities T = M/L, 
TZ = Rh/Re, and <S = ah/cr. Focusing first on the 
edge-on FP, from equations (4) and (7) we obtain: 

T 

^ „ p0.2n.42 /o\ 

7^1.1 51.73 °^ ■ («J 

Thus, in order to transform the DM halos FP into 
the observed FP, the product T7?.^^-^<S~^'^^ must 
systematically increase as R^-^a^-"^^, which in turn, 
again from equation (4), is approximately propor- 
tional to In principle, T, TZ, and S could all 
vary in a combined and regular way from cluster to 
cluster, so that equation (8) is satisfied. Of course, 
given the small scatter around the best fit rela- 
tion (4), this kind of solution requires a remarkable 
fine tuning of the variations of the three parame- 
ters. Alternatively, it is possible that only one of 
the three parameters varies significantly, while the 
other two are approximately constant. This situa- 
tion is analogous to that faced in the studies of the 
physical origin of the FP tilt of elliptical galaxies, 
where the so called "orthogonal exploration of the 
parameter space" is often adopted (see, e.g., Ren- 
zini & Ciotti 1993; Ciotti 1997) . In this approach 
all but one of the available model parameters are 
fixed to constant values, and the goal is to deter- 
mine what kind of variation of the "free" parame- 
ter is necessary to reproduce the observed FP tilt. 
Several quantitative results have been derived in 
this framework (see, e.g., Ciotti, Lanzoni & Ren- 
zini 1996; Ciotti & Lanzoni 1997, and references 



therein), even though, by construction, it cannot 
provide the most general solution to the problem, 
and the choice of the specific parameter responsi- 
ble for the tilt is somewhat arbitrary (see Bcrtin ct 
al. 2002; Lanzoni & Ciotti 2003). In the present 
context some of this arbitrariness can be removed: 
in fact, here we assume that 1) the DM distribu- 
tion in real clusters is described by the simulated 
DM halos, and that galaxies are merely dynamical 
tracers of the total potential well, 2) in addition to 
the edge-on FP, we also consider the constraints 
imposed by the FJ and the Kormendy relations. 
These two points will allow to use the orthogonal 
exploration approach for determining what is the 
most plausible origin of the tilt between the simu- 
lated and the observed cluster FP. 

In order to make the DM halos FP reproduce 
the observed one within the framework of the 
orthogonal exploration approach, we have three 
different possibilities, each corresponding to the 
choice of T, TZ, or <S as the key parameter, while 
keeping constant the remaining two in the l.h.s. 
of equation (8). Note that the two choices based 
on variations of 72. or <S should be interpreted from 
an astrophysical point of view as systematic differ- 
ences in the way galaxies populate the cluster DM 
potential well as a function of the cluster mass. 
However, the orthogonal analysis of the FJ and the 
Kormendy relations strongly argue against these 
two solutions, since from equations (1), (2), (5), 
and (6) one obtains 



and 

T^'^Re''- (10) 

Thus, it is apparent that any attempt to repro- 
duce equation (8) by a variation of TZ (or S) alone 
will fail at reproducing the F J (or the Kormendy) 
relation: in fact, the only parameter appearing in 
all the equations (8), (9), and (10) is the mass-to- 
light ratio^ T. 

Therefore, while a purely structural (TZ) and a 
purely dynamical {S) origin of the tilt between the 

^Note that the constraints imposed by the FJ and the Kor- 
mendy relations should not be considered redundant with 
respect to those imposed by the edge-on FP: in fact, these 
two relations, albeit with a large scatter, describe how 
galaxies are distributed on the face-on FP. 
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Table 3 

Exponents and rms scatter of the scaling relations M oc aj^ , M oc , and M oc i?^^ a' 



fr 


a 


rms 


b 


rms 


c 


d 


rms 


0.8 


3.01 


0.12 


2.67 


0.10 


1.47 


1.43 


0.04 


0.5 


3.02 


0.12 


2.36 


0.15 


1.10 


1.73 


0.04 


0.3 


3.10 


0.13 


2.29 


0.17 


0.97 


1.95 


0.04 


0.1 


3.12 


0.13 


2.15 


0.19 


0.88 


2.02 


0.03 



Note. — M is halo total mass, Rfr is the projected radius 
containing the mass fr X M, and afr is the mean projected 
velocity dispersion within Rfr'. Ho. 5 = Rh- 



DM halos FP and the clusters FP seem to be both 
ruled out by the reasons above, a systematically 
varying mass-to-light ratio, for TZ and S constant, 
could in principle account for all the three consid- 
ered scaling relations. In particular, from equa- 
tion (8), T cx with a ^ 0.3. Guided by this 
indication, we tried to superimpose the points cor- 
responding to the simulated DM halos to the sam- 
ple of observed clusters by using T oc M^, and we 
found that if 



280/1 



M 



1014 Mq/H 



0.23 



^, (11) 



the edge-on FP of DM halos is practically indis- 
tinguishable from that of real clusters (sec Fig. 2a 
and Fig. 3): the value of a derived from this as- 
sumption is a — (3 /{I — (3) ~ 0.3. Note that such 
a trend of T is in agreement not only with the; 
expectations of equation (8), but also with what 
inferred observationally in the B-band (S93; Gi- 
rardi et al. 2002; Bahcall & Comeford 2002; but 
see Bahcall, Lubin & Dorman 1995; Bahcall et 
al. 2000; and Kochanek et al. 2003 for claims 
of a constant mass-to-light ratio at large scales), 
and with what inferred from the comparison be- 
tween the observed B-band luminosity function of 
virialized systems and the halo mass function pre- 
dicted in CDM cosmogonies (Marinoni & Hudson 
2002). A remarkable agreement is also found with 
the results of van den Bosch, Yang & Mo (2003), 
who, using a completely different approach, ob- 
tain T a M"'^^ for their model A, which al- 
lows for a non-constant T at the cluster scales 
(see their equation 15 and Table 2). In addition, 
T ~ 280/iMq/Lq for M ~ IO^^Mq/Zi clusters, 
and T ~ 475/iMq/L0 for M ~ IO^^Mq/Zi clus- 



ters are values well within the range of the various 
estimates for galaxy clusters (e.g., Adami et al. 
1998b; Mclhcr 1999; Wilson, Kaiser & Luppino 
2001; Girardi et al. 2000; Girardi et al. 2002). 
It is also remarkable (as not a necessary conse- 
quence) that by adopting equation (11) also the 
face-on FP, and the FJ and Kormendy relations 
are very well reproduced, as apparent from Fig. 2b 
and Figs. lab, respectively. 

What could be the physical interpretation of 
the required trend of the mass-to-light ratio with 
cluster mass? Wc note that equation (11) can 
be formally rewritten as Tgai x (M/Mgai) oc M^, 
where M is the DM mass of the clusters (^ their 
total mass), Mgai is their total stellar content in 
galaxies. 



^ / jVgai(£)T,(L)LdL 



(12) 



where Ngai{L) is the cluster luminosity function, 
and finally T*(L) is the mean stellar mass-to-light 
ratio of a galaxy of total luminosity L. Thus, the 
trend of T with M could be ascribed or to Tgai oc 
Ml^ for M/Mgai = const, or to Mgai oc M^-^^ 
for Tgai = const (or, more generally, to a com- 
bined effect of these two quantities). Both pos- 
sible solutions have interesting astrophysical im- 
plications. For example, a constant Tgai from 
cluster to cluster is obtained only if all clusters 
have the same population of galax;ies, i.e., if they 
are characterized by a universal luminosity func- 
tion (LF) and by a similar morphological mix, so 
that the distribution of the stellar mass-to-light 
ratios of their galaxies is also the same. In such a 
case, the required trend Mgai oc M^-"^"^ should be 
entirely explained by a systematic increase with 
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M of the total mimbcr of galaxies, in remarkable 
agreement with several studies of the halo occupa- 
tion numbers, that find -/Vgai oc M", with a in the 
range 0.7 0.9 (Peacock & Smith 2000; Scranton 
2002; Berhnd & Weinberg 2002; Marinon & Hud- 
son 2002; van den Bosch et al. 2003). To study the 
ratio M/Afgai, an estimate of the stellar mass in 
galaxies can also be obtained from the observed to- 
tal luminosity in the near infrared (since Es are the 
dominant component of the cluster population and 
their luminosity in the K band gives a reasonable 
measure of their stellar mass, i.e., Mgai oc Lk)- 
However, observational results arc still uncertain 
and controversial: a constant or weakly decreasing 
M/Lk with M is reported by Kochaneck et al. 
(2003) for the 2MASS clusters, while an increas- 
ing M/Lk is claimed by Lin, Mohr & Standford 
(2003) for the data from the same survey. 

In any case, the universality of the cluster lu- 
minosity function is still under debate since it ap- 
pears to be appropriate in many cases, but vari- 
ations in some individual clusters have also been 
reported (see, e.g., Yagi et al. 2002, De Propris 
et al. 2003, and Christlein & Zabludoff 2003 for 
detailed discussions and recent results). In partic- 
ular, the LF of early-type (or quiescent) galaxies 
appears to vary significantly among clusters (but 
see Andrcon 1998): if the richer clusters contain a 
proportionally larger fraction of elliptical galaxies 
(Balogh et al. 2002), that are characterized by a 
higher T* compared to that of spirals, a cluster- 
dependent Tgai should be expected even in pres- 
ence of a universal LF. 

5. Cluster vs. galaxy scaling relations 

As discussed in Section 2, the cosmological col- 
lapse model predicts the same scaling relations at 
all mass scales and, remarkably, the edge-on FP 
and the (although more dispersed) Kormendy re- 
lation of clusters and ellipticals are very similar. 
In particular, the interpretation of the FP tilt of 
elliptical galaxies as the combined results of the 
virial theorem and strong structural and dynam- 
ical homology, implies T* cx L^ "^ in the B band 
(e.g., Fabcr et al. 1987; van Albada, Berlin & 
Stiavelli 1995; Berlin et al. 2002), in agreement 
with the result for clusters (equation 11). In Sec- 
tion 4 we showed that equation (5) coupled with 
T cx L^ "^ does reproduce the observed FJ relation 



at cluster scales, but clearly this cannot work in 
the case of galaxies: in fact, the FJ relation of 
Es is characterized by a significantly higher expo- 
nent (~ 4) than that of clusters (^ 2). Indeed, 
T* should actually decrease for increasing galaxy 
luminosity in order to reproduce the observed FJ, 
at variance with all the available indications (e.g., 
Faber et al. 1987; J0rgensen et al. 1996; Scodeg- 
gio et al. 1998). 

One is then forced to assume that 1) the rela- 
tion M oc cr^ij, does not apply in the case of Es, per- 
haps due a failure of standard cosmology at small 
scales, or 2) evolutionary processes have modified 
the ratio cr*/avir, where cr* is the galaxy central 
velocity dispersion (from which the FJ relation is 
derived). For what concerns point 1), this might 
be another problem encountered by the current 
cosmological paradigm at small scales, in addition 
to the well known over-prediction of the number of 
DM satellites (the so called "DM crisis" ; see Moore 
et al. 1999; but see also Stohr et al. 2002 and 
the recent results on the "running spectral index" 
obtained by Peiris et al. 2003 from the WMAP 
data) and the cusp-core problem for low surface 
brightness and dwarf galaxies (e.g., Swaters et al. 
2003ab, and references therein). However, this 
point is at present very speculative, and thus in 
the following we restrict to the standard CDM 
paradigm, that predicts M oc a^^^ also at galac- 
tic mass scales. Two physical processes (namely, 
gas dissipation and early-time merging) certainly 
played a major role in galaxy evolution, thus sup- 
porting point 2). The effects of these two processes 
have been already discussed in the context of k- 
space by Burstein et al. (1997): here we will focus 
on the implications that can be derived from the 
simpler FJ relation, taking into account the addi- 
tional information from recent numerical simula- 
tions (NLC03ab) and the constraints imposed on 
galaxy merging by the recently discovered Mbh-c 
relation (Ferrarese & Merritt 2000; Gebhardt et 
al. 2000), according to which Mbh oc c7~^, with 
very small dispersion. 

We then adopt the point of view that at the 
epoch of the detachment from the Hubble flow, 
also the seed galaxies were characterized by the 
"universal" scaling law M oc a^j^., and we discuss 
the possibility that galaxy merging and gas dissi- 
pation originated the observed FJ of Es. We start 
by considering the effect of dissipationless merg- 
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ing alone. From this point of view, it is clear that 
merging played different roles in the evolutionary 
history of clusters and galaxies. In fact, while clus- 
ters can bo thought as formed from the collapse of 
density perturbations at scales that just became 
non-linear, galaxies are presently in a highly non- 
linear regime and the merging they suffered since 
their separation from the Hubble expansion can no 
longer be interpreted in terms of the cosmological 
collapse of density fluctuations. The differences 
between these two dynamical processes have im- 
portant consequences for the present discussion. 
In fact, in the cosmological collapse case the ini- 
tial conditions correspond to those of a "cold" sys- 
tem (i.e., 2Ti + Ui = V < 0) , and so virialization 
increases the virial velocity dispersion of the end- 
products as Tf =Ti — V: the systematic increase 
of CTvir with M in clusters is basically due to this 
process. At highly non-linear scales (such as, for 
instance, in the case of galaxies in the outskirts 
of clusters or groups) the situation is considerably 
different. In fact, if merging occurs, it cannot be 
interpreted as the collapse of a cold system, but, 
on the contrary, the initial conditions of the merg- 
ing pair arc in general characterized by a null or 
a positive V (see, e.g., Binney & Tremaine 1987, 
Chap. 7): under these conditions, the virial veloc- 
ity dispersion of the remnant will not increase. 
In fact, numerical simulations of one and two- 
component galaxy models show that successive 
dissipationless merging at galactic scales, while 
preserving the edge-on FP, does not reproduce the 
FJ relation, being the end-products characterized 
by a too low (i.e.. nearly constant) ct, compared to 
the expectations of the FJ (NLCOSab). In other 
words, dissipationless merging at galactic scales 
in general will produce a relation M oc a"^^ with 
a > 3. Obviously, this could be an interesting 
property in our context, since we are exactly look- 
ing for a mechanism able to increase a from 3 to 
4. We also note that faint elliptical galaxies do 
follow a FJ relation with a best-fit slope signifi- 
cantly lower than 4 (Davies ct al. 1983 report a 
value of ~ 2.4) and more similar to the cosmo- 
logical predictions: a simple interpretation of this 
fact could be that small galaxies experimented less 
merging events than giant Es, so that their scaling 
relations are more reminiscent of the cosmological 
origin. However, even if able to increase the expo- 
nent of the M-a relation, several theoretical and 



empirical arguments clearly indicate that purely 
dissipationless merging cannot be at the origin of 
the spheroids. For example, why the exponent 
should increase to the observed value is unclear, 
even though it has been claimed that the F J of Es 
is the result of cumulative effects of inelastic merg- 
ing and passing encounters taking place in cluster 
environment (Funato, Makino & Ebisuzaki 1993; 
Funato & Makino 1999). These authors already 
discuss some weakness of the scenario they sug- 
gest (for instance, the difficulty to account for the 
analogous Tully-Fisher relation of spirals), and we 
add that it is not clear whether also the other scal- 
ing laws are reproduced by their simulations, nor 
how the FJ relation could be followed also by the 
field Es if resulting from dynamical interactions 
in cluster environment. Moreover, a further diffi- 
culty is added by the recently discovered Mbh-ct 
relation: if the FJ is the cumulative result of ran- 
dom dynamical processes that changed the galaxy 
velocity dispersion, how can the Mbh-o" relation 
be so tight? Finally, an even stronger empirical 
argument argues against substantial dissipation- 
less merging for Es: in fact, if the merging end- 
products arc forced to obey both the edge-on FP 
and the Mbh-ct relation, then they are character- 
ized by exceedingly large effective radii with re- 
spect to real galaxies (Ciotti & van Albada 2001). 

Concerning the effect of gas dissipation on the 
predicted relations at galactic and cluster scales, 
we note that an important difference already re- 
sides in what is effectively observed in the two cases 
when deriving the scaling laws: the galaxy distri- 
bution, tracing the gravitational potential of the 
dark matter, in the case of clusters; the stellar 
population within a fraction of the effective ra- 
dius, where stars themselves arc the major contrib- 
utors to the total gravitational field, in the case of 
Es. Since the top-hat model properly describes the 
(dissipationless) collapse of the DM component, it 
is not surprising that it cannot be accurately ap- 
plied for predicting the properties of the baryonic 
(dissipative) matter where it is dominant. In ad- 
dition, gas dissipation is thought to be negligible 
for the formation and evolution of clusters, while 
it empirically appears to be increasingly impor- 
tant with mass for the galaxies, as mirrored by 
their Mg2-a and color-magnitude relations, and 
by the observed metallicity gradients (e.g., Saglia 
et al. 2000; Bernardi et al. 1998; Bernardi et 
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al. 2003a). As discussed by Ciotti & van Albada 
(2001), its effects also represent a possible solu- 
tion to the exceedingly large effective radii found 
for (dissipationlcss) merger products forced to fol- 
low both the FP and the Mbh-c relation. Thus, 
from the considerations above, gas dissipation is a 
good candidate to explain the difference between 
cluster and galaxy FJs. However, we note that 
baryonic dissipation alone cannot solve the dis- 
crepancy: in fact, its effect is to increase ci*, thus 
further decreasing the exponent in the galactic FJ 
relation. This can be seen in a more quantitative 
way as follows. Given that galactic DM halos fol- 
lows M oc a^jj., then the FJ relation implies that 
(M/M*)T* ((T*/crvir)^ OC a~^: this means that at 
fixed M/M^, and T*, dissipation should have been 
more important in low mass galaxies, contrarily to 
the existing observational evidences. 

In summary, the arguments presented above 
strongly suggest that merging and gas dissipation 
played a fundamental role in the formation and 
evolution of galaxies. In particular, given the com- 
petitive effect of these two processes on the slope of 
the M-cr relation, they appear to be both necessary 
for modifying the M oc ct^jj. into the observed FJ. 
Coupled with the necessity of dissipative merging, 
the observational evidence that the bulk of stars in 
Es, even at high z, arc old puts a strong constraint 
on the time when the mass of spheroids was assem- 
bled: say, z>2. As already recognized by several 
authors, this is in agreement with the available in- 
formation about the star formation history of the 
Universe and the redshift evolution of the quasar 
luminosity function (see, e.g., Haehnelt & Kauff- 
mann 2000; Ciotti & van Albada 2001; Burkert & 
Silk 2001; Yu fc Tremaine 2002; Cavaliere & Vit- 
torini 2002; Haiman, Ciotti & Ostriker 2003). A 
tight connection between the star formation and 
the quasar activity in Es is further supported by 
the striking similarity of slopes between the FJ 
and the Mbh-c relations, thus suggesting a com- 
mon origin for these two empirical laws. 

Based on these considerations, it seems likely 
that the scaling relations of Es are the result of the 
cosmological collapse of the density fluctuations 
at the epoch when galactic scales became non- 
linear, plus important modifications afterward due 
to early-time merging in gas-rich systems. On the 
contrary, since the present day non-linearity scale 
is that of galaxy clusters, and since the role of gas 



dissipation is thought to be marginal in clusters 
life, the cosmological collapse model and a mass- 
to-light ratio increasing with cluster mass seem to 
be sufficient to account for their scaling relations. 

6. Discussions and conclusions 

In this paper we have addressed the question 
of whether high mass DM halos (that are thought 
to host the clusters of galaxies) do define the FP 
and the other scaling relations observed for nearby 
clusters. For that purpose, we have analyzed a 
sample of 13 massive DM halos, obtained with 
high-resolution N-body simulations in a ACDM 
cosmology. After verifying that the DM halos 
do follow the predictions of the spherical collapse 
model for virialized systems, we have found that 
also their projected properties define a FJ, a Ko- 
rmendy and a FP-like relations. This latter re- 
sult is not trivial and it can be traced back to the 
weak (structural and dynamical) homology shown 
by the DM halos assembled via hierarchical cosmo- 
logical merging. However, the slopes of the DM 
halos scaling laws do not coincide with the ob- 
served ones, and we have discussed what kind of 
systematic variations of one or more of the struc- 
tural and dynamical properties of the galaxy dis- 
tribution with cluster mass could account for such 
differences. We have shown that two of the three 
basic options can be discarded just by requiring 
the simultaneous reproduction of the (edge-on) 
FP, FJ and Kormendy relations. A solution that 
instead works remarkably well is to assume a clus- 
ter mass-to-light ratio T increasing as a power law 
of the luminosity. The required normalization and 
slope well agree with those estimated observation- 
ally for real galaxy clusters. We have discussed 
two possible causes for such a trend of T, namely 
a systematic increase of the galaxy mass-to-light 
ratio with cluster luminosity, or a decrease of the 
baryonic mass fraction for increasing cluster total 
mass, concluding that both the possibilities are 
consistent with the available observational data. 
In any case, it appears that the FJ, Kormendy and 
FP relations of nearby clusters of galaxies can be 
explained as the result of the cosm,ological collapse 
of density fluctuations at the appropriate scales, 
plus a systematic trend of the total mass-to-light 
ratio with the cluster mass. Note that this is by no 
means a trivial result: in fact, it is known that nu- 
merical simulations of galaxy dissipationless merg- 
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ing (where the mass-to-light ratio is kept fixed) do 
reproduce well the edge-on FP, but badly fail with 
the FJ and the Kormendy relations (NLC03b). 

We next focused on the fact that Es follow a FJ 
relation with a significantly different slope (~ 4) 
compared to that of clusters 2), but show a 
similar trend of the mass-to-light ratio with the 
system luminosity. We have discussed the impli- 
cations of this observational evidences under the 
assumption that at the epoch of the; detachment 
from the Hubble flow, M oc cr^jj, also at galac- 
tic scales, as required by standard cosmology, and 
taking into account several empirical and theo- 
retical evidences suggesting that gas dissipation 
and merging must have played an important role 
in galaxy evolution. Since dissipationless merging 
and baryonic dissipation have competitive conse- 
quences on the system velocity dispersion, we ar- 
gue that combined effects of these; two processes 
are required to account for the slope of the FJ of 
elliptical galaxies. We therefore conclude that the 
scaling relations of Es might be the result of the 
cosmological collapse of density fluctuations at the 
epoch when galactic scales became non-linear, plus 
successive modifications due to (early-time) dissi- 
pative merging. This scenario seems to be sup- 
ported also by the similarity between the slopes 
of the FJ and the Mbh-c relations, and between 
the peak of the star formation rate of the Universe 
and that of the quasar luminosity function. 

Before to conclude we point out that, while no 
observational data are yet available for the scaling 
relations of galaxy clusters at high redshift, the 
results that we have obtained allow us to specu- 
late on what could be expected. On one hand, 
the predictions of the top-hat model arc the same 
at any redshift, and also weak homology of DM 
halos in numerical simulations appears to hold at 
any epoch: thus, if the dependence of the cluster 
T on i changes accordingly to passive evolution, 
scaling relations with the same slopes should be 
found also at high redshift. Even if with some 
uncertainties, observational evidences for passive 
evolution of these cluster properties already exist. 
In fact, the characteristic luminosity of the cluster 
LF increases with redshift consistently with pure, 
passive luminosity evolution models (De Propris 
et al. 1999), and the same is found for the stel- 
lar mass-to-light ratio of elliptical galaxies, as 
derived from the studies of their FP up to red- 



shift z = 1.27 (e.g., van Dokkum & Franx 1996; 
J0rgensen et al 1999; Kelson et al. 2000; Treu, 
M0ller & Bertin 2002; van Dokkum & Standford 
2003). On the other hand, a more complicate evo- 
lution of T = Tgai X M/Mgai with z might also be 
expected, and thus the redshift variation of the FP 
zero-point and slope could give important infor- 
mation on cluster and galaxy evolution. Although 
the fraction of Es in clusters appears not to change 
with redshift, a substantial increase of the relative 
amount of spirals with respect to SO galaxies at 
higher z is reported (Dressier et al. 1997; Fasano 
et al. 2000). In addition, also M/Mgai might in- 
crease with z, if clusters form by continuous accre- 
tion of galaxies from the field, in a DM potential 
well that already settled on shorter time scales. 
Of course, a detailed prediction of the evolution of 
the cluster scaling relations in this latter case is 
much more difficiilt. 

Due to the interest of the questions addressed 
above, it is apparent that with more available 
data (e.g., from the SDSS, RCS, MUNICS, 2dF, 
2MASS surveys), a strong effort should be devoted 
to better determine the cluster scaling relations 
while, on theoretical ground, a larger set of sim- 
ulated clusters, spanning a wider mass range and 
taking into account also the presence of gas, should 
be performed and analyzed for different choices of 
the cosmological parameters. 
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log a log 

Fig. 1. — Panel a: FJ relation for the observed clusters (filled triangles) and the DM halos when equation 
(11) is used for the mass-to-light ratio (empty circles). The corresponding best fit relations are L oc a"^-^^ 
(dashed line) and L cx cr^'^^ (solid line), respectively. Panel b: The Kormcndy relation for the same data as 
in Panel a. The best fit relations are L oc Rl'^^ for the data (dashed line), and L oc Rl'^^ for the DM halos 
(solid line). Each DM halo is represented by 3 empty circles corresponding to the 3 line-of-sight directions. 
Luminosities are normalized to IO^^Lq/Zi, velocity dispersions are in 1000 km/s, and radii in Mpc/h. 
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Fig. 2. — The edge-on (upper panel) and facc-on (lower panel) views of the FP for the observed clusters (filled 
triangles) and the DM halos when equation (11) is used for the mass-to-light ratio (empty circles). The solid 
line would be the best fit of the DM halos when assuming a constant mass-to-light ratio T = 28O/1M0/L0. 
Note that mass increases for decreasing pi . 
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Fig. 3. — The FP of observed clusters (filled triangles) and DM halos when equation (11) is used for 
mass-to-light ratio (empty circles). 
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